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Abstract 

The  establishment  of  flow  conditions  character¬ 
istic  of  hypersonic  flight  in  a  ground  test  facility  by 
conventional  means  requires  plenum  conditions 
that  have  very  high  temperatures  and  give  rise  to  a 
variety  of  materials  problems.  A  well-known 
approach  to  solving  these  problems  Is  to  use  high- 
pressure,  moderate  temperature  plenum  condi¬ 
tions  and  add  energy  as  the  flow  expands  through 
a  carefully  designed  nozzle.  A  modern  approach 
(Radiatively  Heated  Wind  Tunnel)  uses  an  HF 
laser  to  excite  the  (02^1)  mode  of  CO2  naturally 
present  or  added  to  the  flow.  One  of  the  early  chal¬ 
lenges  in  the  development  will  be  to  demonstrate 
the  ability  to  add  energy  into  the  supersonic  flow  in 
a  predictable  and  controllable  manner.  Conse¬ 
quently,  a  subscale  experiment,  the  Laser  Demon¬ 
stration  Device  (LDD),  is  planned  to  investigate  the 
transfer  of  laser  energy  to  a  small,  expanding  noz¬ 
zle  flow. 

Previous  modeling  of  the  energy  transfer  of  the 
laser  energy  into  the  gas  for  the  LDD  depended  on 
a  conventional  two-temperature  approach  to 
describe  the  internal  structure  of  the  molecules  and 
to  determine  the  rate  coefficients.  Realizing  that 
the  detailed  physics  of  the  energy  transfer  was 
important  to  the  description  of  the  wind  tunnel  per¬ 
formance  and,  hence,  its  viability,  an  effort  was 
undertaken  to  compute  the  effect  of  CO 2  relaxation 
using  a  technique  in  which  the  individual  vibrational 
states  are  considered.  Results  of  these  computa¬ 
tions  are  reported  here. 

The  state-to-state  approach  has  been  utilized 
before  in  rocket  exhaust  flows  where  the  densities 
are  low  and  temperatures  modest.  Individual  states 


of  CO2,  N2,  and  O2  were  maintained  as  distinct 
states,  and  relaxation  rates  for  each  state  were 
determined.  A  vibrational  temperature  is  not 
required.  The  laser  excitation  along  the  flow  direc¬ 
tion  was  described  by  an  empirical  absorption 
function  and  Is  the  same  as  used  previously.  Feed¬ 
back  coupling  between  the  laser  absorption  func¬ 
tion  and  the  state  distributions  is  not  included  at  the 
present  level  of  approximation. 

The  results  shown  here  are  in  agreement  with 
earlier  work;  significant  energy  coupling  to  the  flow 
and  laser  heating  of  the  flow  appears  to  be  a  feasi¬ 
ble  approach.  However,  the  details  of  the  vibra¬ 
tional  relaxation  show  that  a  conventional  two-tem¬ 
perature  approach  does  not  describe  the  resultant 
distribution.  Consequently,  inclusion  of  these 
effects  will  be  required  for  successful  prediction  of 
exit  plane  properties,  and  thus,  will  be  of  critical 
importance  to  facility  design  and  control. 

Introduction 

The  development  of  hypersonic  flight  vehicles 
Is  seriously  limited  by  the  lack  of  ground  test  facili¬ 
ties  that  can  routinely  produce  flows  in  excess  of 
Mach  8.  At  these  higher  flight  Mach  numbers,  res¬ 
ervoir  conditions  in  conventional  expansion-type 
wind  tunnels  become  extreme  enough  that  con¬ 
tainment  of  the  reservoir  becomes  a  limiting  factor 
on  the  tunnel  operation.  For  example,  to  achieve  a 
simulated  flight  condition  of  Mach  10  at  30-km  sim¬ 
ulated  altitude  (0.03  atm,  200  K),  a  simple  calcula¬ 
tion  shows  that  chamber  pressure  in  excess  of 
1,000  atm  and  temperatures  in  excess  of  4,000  K 
are  required.  This  temperature  is  well  above  the 
melting  temperature  of  materials  that  can  be  used 
for  wind  tunnel  construction.  Additionally,  at  these 
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high  temperatures  there  is  significant  air  chemistry, 
and  the  working  gas  delivered  to  the  test  article  Is 
not  representative  of  the  air  at  actual  flight  condi¬ 
tions.  Further,  at  the  high  pressures,  throat  erosion 
can  significantly  degrade  flow  quality,  even  at 
lower  temperatures.  One  approach  to  achieving 
these  high  Mach  number  flows  has  been  to  use 
facilities  that  operate  for  very  short  times  so  that 
materials  degradation  is  minimal  and  accept  what¬ 
ever  chemical  effects  are  present.  This  approach 
limits  test  time  for  component  testing  and  the 
chemistry  effects  seriously  limit  the  utility  for  com¬ 
bustor  development  and  testing. 

An  approach  to  attack  this  problem  has  been 
proposed^  in  which  the  energy  Is  added  with  a 
high-power  laser  to  the  supersonic  flow  down¬ 
stream  of  the  throat.  With  this  approach,  the  ple¬ 
num  can  be  at  very  high  pressure  but  at  a  relatively 
cool  temperature  through  the  throat  expansion, 
thus  avoiding  materials  and  air  chemistry  problems 
that  characterize  conventional  approaches.  This 
energy  addition  is  accomplished  through  a  portion 
of  the  flow  In  which  the  expansion  is  relatively  slow 
(although  maintained  at  supersonic  speeds)  and 
then  a  final  expansion  is  accomplished  to  the 
desired  Mach  number  and  test  conditions. 
Although  there  are  other  candidate  methods  for 
energy  addition  (e.g.,  electron  beam,  microwave, 
and  magnetohydrodynamics),  the  use  of  high- 
power  lasers  has  progressed  farthest  in  the  design 
studies  and  is  the  subject  of  this  paper. 

As  described  in  Refs.  1-4,  development  of  such 
a  facility,  the  radiatively  heated  wind  tunnel 
(RDHWT),  presents  several  significant  engineering 
challenges.  Consequently,  incremental  proof-of- 
concept  demonstrations  are  planned  as  part  of  the 
path  toward  development  of  a  full-scale  device.  An 
early  challenge  that  is  critical  to  validating  the  feasi¬ 
bility  of  the  approach  is  the  demonstration  of 
energy  addition  to  a  supersonic  flow  for  an  increase 
In  the  test  chamber  temperature  and  pressure 
under  controllable  and  predictable  conditions.  Such 
a  demonstration  Is  quite  feasible  with  HF  or  DF 
lasers  that  currently  exist,  and  plans  are  proceed¬ 
ing  to  accomplish  just  such  a  demonstration.  One 
approach  for  the  demonstration  is  to  use  a  1 0-kW 
HF  laser,  the  Laser  Demonstration  Device  (LDD), 
with  C02-seeded  airflow.  The  high-power  HF  laser 


already  exists,  and  it  has  been  shown  earlier^  that 
there  are  several  CO2  absorption  lines  that  are 
close  to  the  laser  lines.  In  the  high-pressure  envi¬ 
ronment  planned  for  the  absorption  region,  there 
will  be  significant  coupling  of  the  input  laser  energy 
Into  the  CO2.  Subsequent  transfer  of  the  energy 
from  the  CO2  into  the  kinetic  modes  of  the  flow  was 
shown  to  be  plausible  based  on  a  straightforward 
two-temperature  relaxation  model  in  which  the  N2  - 
O2  V-V  transfer  is  the  limiting  rate  for  the  relaxation. 
Although  the  feasibility  model^  indicated  that  there 
is  significant  coupling  of  the  laser  energy  into  the 
flow,  in  actuality  the  laser  couples  to  specific  vibra¬ 
tional  states  in  the  CO2  and  the  details  of  the  sub¬ 
sequent  relaxation  are  not  modeled  in  the  conven¬ 
tional  approach.  Consequently,  it  was  decided  to 
investigate  the  thermalization  of  the  Input  energy  of 
the  LDD  experiment  using  a  state-specific  relax¬ 
ation  (SSR)  model  that  includes  the  individual 
molecular  vibrational  states.®"® 

The  SSR  model  has  the  capability  to  model  the 
individual  molecular  vibrational  states  participating 
in  the  laser  absorption  and  subsequent  thermaliza¬ 
tion  as  separate  entities.  Thus,  each  vibrational 
state  considered  becomes  a  separate  species.  Its 
concentration  Is  a  result  of  the  kinetic  processes  to 
which  the  molecules  are  subjected,  and  a  state 
specific  NLTE  distribution  is  calculated.  Vibrational 
temperature  is  not  used  explicitly,  but  rather  may 
be  a  consequence  if  the  excited  states  exhibit  the 
properties  characteristic  of  such  a  temperature, 
i.e.,  a  Boltzmann  distribution.  The  SSR  model  has 
been  used  previously  to  investigate  the  details  of 
the  CO2  radiative  transfer  from  rocket  exhausts 
and  Includes  the  vibrational  states  participating  in 
the  present  laser  absorption  and  subsequent  ther¬ 
malization.  An  earlier  study®  showed  that  different 
distribution  models  can  result  in  different  radiative 
transfer  properties,  and  it  was  surmised  that  includ¬ 
ing  the  details  of  the  relaxation  processes  would 
show  some  differences  in  the  predicted  energy 
transfer  to  the  medium. 

In  the  present  paper,  the  elements  of  the  vari¬ 
ous  models  for  modeling  the  laser  addition  to  the 
flow  are  described.  These  elements  include  the 
specification  of  the  LDD  nozzle  and  power  addition 
profile,  the  flow  solver,  the  molecular  model,  the 
rate  coefficients,  and  the  thermodynamics.  Results 
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of  application  to  the  LDD  are  described  and  the 
effect  of  the  energy  addition  on  the  exit  plane  gas 
dynamic  properties  are  examined.  The  internal 
state  distributions  are  examined  for  their  relevance 
to  a  conventional  two-temperature  approach  and 
the  effectiveness  of  the  thermalization.  Implications 
for  the  larger  scale  devices  are  drawn. 

Models 

Laser  Demonstration  Device  Model 

The  nozzle  contour  used  for  the  present  study 
is  shown  in  Fig.  1  and  is  similar  to  a  contour  pro¬ 
posed  early  in  the  effort.  Although  the  exact  con¬ 
tour  of  the  nozzle  is  still  in  the  design  and  develop¬ 
ment  phase,  the  larger  issues  of  examining  the 
thermalization  of  laser  energy  can  be  studied  with¬ 
out  the  exact  contour.  Common  elements  of  all 
nozzle  designs  for  subscale  feasibility  demonstra¬ 
tion  include  the  rapid  expansion  from  the  very 
small  throat  (0.28-cm  diam  here)  to  an 
area  ratio  of  about  2,  followed  by  a  long 
region  of  constant  or  slowly  increasing 
area  ratio  and  then  fairly  rapid  expan¬ 
sion  to  a  larger  area  ratio  (here  about 
90,  exit  diameter  of  about  2.6  cm).  The 
reservoir  conditions  used  for  this  dem¬ 
onstrator  are  80  atm  and  1,150  K.  The 
working  gas  is  air  seeded  with  CO2  to 
bring  the  mixture  to  40-percent  CO2  and 
60-percent  air. 


The  energy  in  the  LDD  is  supplied  with  a  10-kW 
HF  chemical  laser  and,  except  for  power  level,  is 
expected  to  be  similar  to  that  described  for  larger 
devices.^  The  energy  addition  profile  derived  in  the 
same  manner  as  in  Ref.  3  is  shown  in  Fig.  2.  It 
should  be  noted  that  this  is  a  specified  profile  and 
is  an  input  condition  for  the  calcuiations  reported 
here.  Coupiing  between  the  present  detailed  model 
and  the  input  laser  beam  has  not  been  accom¬ 
plished  yet.  As  pointed  out  above,  there  are  sev¬ 
eral  laser  lines  nearly  coincident  with  absorption 
lines  in  the  2.7-pm  band  of  C02-  In  the  LDD  exam¬ 
ined  here,  the  pressure  in  the  energy  addition 
region  is  about  6  -  7  atm.  Thus,  there  will  be  signif¬ 
icant  pressure  broadening  of  the  absorption  lines, 
and  energy  transfer  from  the  laser  beam  to  the 
vibrational  modes  of  CO2  will  be  effective.  The 
laser  beams  are  surmised  to  have  their  maximum 
transfer  of  energy  to  the  gas  in  the  relatively  high- 
pressure  region  near  the  exit  of  the  slowly  varying 
section  at  an  axial  position  of  about  20  cm.  The 


It  is  recognized  that  this  model  for  an 
LDD  nozzle  with  an  area  ratio  of  90  gives 
an  expansion  to  a  Mach  number  of  only 
6  with  no  energy  addition,  which  is  far 
from  what  is  required  for  a  hypersonic  facility. 
However,  the  near-term  objectives  are  for  the 
demonstration  of  energy  addition  into  a  super¬ 
sonic  stream  in  a  predictable  and  controllable 
manner.  Use  of  these  more  modest  conditions 
fits  well  within  the  present  experience  base 
and  will  be  relatively  straightforward  to  accom¬ 
plish.  In  this  fashion,  feasibility  of  the  concept 
can  be  examined  without  requiring  significant 
advances  in  the  state  of  the  art,  and  the  ability 
to  predict  the  exit  plane  conditions  can  be 
tested  as  well  as  developing  the  techniques  for 
controlling  the  test  chamber  conditions. 


0  10  20  30  40 

Axial  Position,  cm 

Fig.  1 .  Nozzle  contour. 
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region  from  about  10  to  15  cm  serves  as  a  high- 
pressure  region  that  absorbs  the  remainder  of  the 
energy  not  absorbed  further  downstream. 

Molecular  Model 

Critical  to  the  understanding  of  the  thermaliza- 
tion  of  the  energy  that  will  be  added  to  the  flow  of 
the  RDHWT  is  the  understanding  of  the  distribution 
of  the  internal  energy  states  of  the  molecules  par¬ 
ticipating  in  the  thermalization.  Nitrogen  and  O2  are 
diatomic  molecules  with  relatively  simple  vibra¬ 
tional  structure,  while  CO2  is  a  triatomic  molecule 
with  a  more  complicated  structure.  It  Is,  however,  a 
linear  molecule  and,  consequently,  there  are  regu¬ 
larities  in  the  term  structure  that  simplify  modeling. 

A  partial  term  diagram  for  CO2  is  shown  in  Fig. 
3.  Only  the  V3  =  0  levels  through  v-|  =  4  and  V2  =  8 
are  shown.  States  from  the  other  series,  V3  =  1  and 
V3  =  2,  are  included  in  the  flow  calculations  and  are 
included  in  the  discussion  below,  but  this  partial 
term  diagram  will  be  used  to  explain  the  complexity 
of  the  CO2  molecule  and  the  simplifying  assump¬ 
tions  used  here.  The  numbers  in  parentheses 
below  each  line  in  the  CO2  portion  of  the  term  dia¬ 
gram  are  the  vi,  V2,  and  V3  (symmetric  stretch, 
bending,  and  asymmetric  stretch)  quantum  num¬ 
bers,  respectively.  The  well- 
known  attribute  of  CO2,  the  "acci¬ 
dental"  degeneracy,  or  Fermi  res¬ 
onance,^  in  the  v-|  and  V2  modes 
(every  second  V2  is  in  near  reso¬ 
nance  with  a  v-j),  is  apparent  in 
Fig.  3.  These  states  (e.g.,  (10^0), 

(02^0)),  are  close  together  ener¬ 
getically,  and  there  is,  in  a  quan¬ 
tum  mechanical  sense,  a  mixing 
of  the  wave  functions  of  these 
states.  Pragmatically,  these 
states  are  virtually  indistinguish¬ 
able  and  it  is  very  easy  to  main¬ 
tain  an  LTE  distribution  of  densi¬ 
ties  among  these  states  that 
exhibit  this  Fermi  resonance. 

The  CO2  states  that  exhibit 
this  resonance  and  close  energy 
coupling,  e.g.,  the  (06®0),  (06^0), 

(1.4^0).  (Oe^O),  (142o),  (22^0), 


(06°0),  (14°0),  (22°0),  and  (30°0)  states  at  about 
4,000  cm“\  are  grouped  together  and  are  treated 
as  a  single  state  for  this  work.  They  are  designated 
with  a  capital  L  in  the  notation,  e.g.,  (06L0).  Those 
states  in  this  grouping  that  are  not  connected 
through  the  Fermi  resonance  are  very  close 
together  energetically  and,  as  will  be  seen  in  a  dis¬ 
cussion  of  the  rate  coefficients,  can  be  expected  to 
exhibit  a  thermal  distribution  among  themselves. 
Grouping  the  states  in  this  manner  reduces  the 
number  of  individual  states  included  on  the  flow 
calculation  to  a  manageable  number  without  sacri¬ 
ficing  the  essential  physics  necessary  to  Investi¬ 
gate  the  details  of  the  thermalization  of  the  input 
laser  power.  There  is  similar  structure  with  each  V3 
value  (V3  =  2  was  the  highest  Included  here),  and 
the  grouping  of  the  v*|  and  V2  states  for  each  V3 
was  accomplished. 

Figure  4  shows  the  resulting  term  diagram  that 
was  used  for  the  calculations  here,  including  the 
N2  and  O2  vibrational  states.  For  the  CO2  modes 
there  are  additional  upper  states  designated  “UO,” 
“U1 ,”  and  “U2.”  Included  in  each  of  these  states  are 
the  states  between  the  highest  individual  state  con¬ 
sidered  and  the  U  state  for  V3  =  0, 1 ,  and  2,  respec¬ 
tively.  These  states  are  connected  colllsionally  to 
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Fig.  3.  CO2  V3  =  0  term  diagram. 
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the  highest  individual  state  for  each  V3  with  a  colli¬ 
sion  number  of  10.  Included  also  are  a  number  of 
states  with  a  “U”  in  their  designation.  These  repre¬ 
sent  manifolds  of  states  above  the  energy  at  which 
individual  states  (or  collection  of  states)  were  kept 
distinct.  Included  In  the  “U”  states  for  each  mole¬ 
cule  are  all  the  states  higher  than  10,000  cm""^  . 
These  states  provide  a  reservoir  of  upper  states 
that  are  connected  to  the  next  lower  states  through 
simple  thermal  collisions  with  a  collision  number  of 
10,000.  These  states  thus  provide  a  reservoir  to 
include  the  rest  of  the  molecular  states  necessary 
for  completeness  and  connect  them  to  the  lower 
states  In  a  manner  that  is  physically  consistent  with 
observations.  With  these  considerations,  it  is  seen 
that  the  modeling  includes  23  distinct  states  of 
CO2,  5  of  N2,  and  5  of  O2. 

It  is  seen  in  Fig.  4  that  the  states  of  N2  are  in 
near  resonance  with  the  (00®V3)  modes  of  CO2, 
whereas  the  excited  states  of  O2  are  at  somewhat 
lower  energies  than  either  the  N2  or  CO2.  The  N2  - 
CO2  and  N2  -  O2  V-V  transitions  are  a  well-known 
pathway  for  the  de-excitation  of  CO2  and  N2.  Inclu¬ 
sion  of  these  states  is  an  important  path  for  the 
thermalization  of  the  input  energy,  since  the  N2  will 
provide  a  path  for  energy  transfer  from  the  CO  2  V3 
excited  modes,  and  the  O2  will  then  provide  a  path 
for  the  de-excitation  of  the  N2. 


The  transition  for  the  laser  absorption  in  also 
shown  in  Fig.  4.  The  proposed  LDD  will  deposit 
energy  into  the  flow  through  several  laser  lines 


C02(U) 


N2(U) 


Fig.  4.  Species  term  diagram. 


between  about  3,500  cm“^  and  3,700  cm“^  These 
lines  are  in  the  spectral  region  for  several  absorp¬ 
tion  lines  from  the  CO2  ground  states  and  the 
(10®1)-(02'l)  states,  which  are  in  the  (02L1)  state. 
Although  the  laser  lines  are  very  narrow  and  a 
Doppler  shift  of  the  absorbing  lines  may  be  signifi¬ 
cant,  absorption  occurs  in  the  higher  pressure  por¬ 
tions  of  the  flow  where  the  rotational  lines  are  pres¬ 
sure  broadened.  It  should  be  noted  that  there  will 
be  an  additional  absorption  transition  from  the 
(001)  state  to  the  (02L2)  state,  the  (002)  to  the 
(02L3)  state,  and  so  on.  Any  final  modeling  of  the 
radiative  heating  will  have  to  include  these  transi¬ 
tions  as  well. 

Flow  Solver 

As  in  the  earlier  work®*®’®""'®  the  One  Dimen¬ 
sional  Kinetic  (ODK)  code**''  Is  used  as  the  flow 
solver  for  this  calculation.  The  ODK  code  calcu¬ 
lates  the  inviscid,  chemically  reacting  nozzle 
expansions  of  gaseous  exhaust  mixtures.  The  spe¬ 
cific  version  used  here  was  included  as  a  subpro¬ 
gram  to  the  larger  Two-Dimensional  Kinetic  Analy¬ 
sis  Program  (TDK)**^  and  allowed  up  to  150  chemi¬ 
cal  reactions  and  40  species.  A  modified  version  of 
the  NASA  Equilibrium  code**®  is  also  included  as  a 
subprogram  of  the  larger  TDK  code.  A  complete 
description  of  the  codes  is  included  in  the  TDK 
manual.^ ^  The  results  of  the  equilibrium  code  for 
the  combustion  chamber  are  used  as  initial  values 
at  a  specified  initial  contraction  ratio.  The  flow  is 
described  by  using  a  defined  pressure  distri¬ 
bution  from  the  combustion  chamber  to  the 
sonic  point  and  geometry-defined  relations 
through  the  supersonic  portion.  The  pres¬ 
sure  distribution  for  the  subsonic  flow  is 
obtained  from  the  equilibrium  solution.  The 
species  production  equations  from  the 
chemistry  are  coupled  to  the  gas-dynamic 
equations  in  the  conventional  way.  The 
resulting  differential  equations  are  inte¬ 
grated  by  a  second-order  implicit  method 
derived  in  Nickerson,  et  al.*'^  Step-size  con¬ 
trol  is  included  to  ensure  that  the  rates  of 
change  of  the  sundry  gas-dynamic,  species 
production,  and  energy  production  terms  are 
maintained  within  the  form  of  the  solution 
-  implied  by  the  second  order  numerical  tech¬ 
nique  used. 
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In  principle,  there  is  no  requirement  that  the 
species  be  distinct  chemical  species.  Treating 
each  excited  vibrational  state  as  a  separate  spe¬ 
cies  allows  one  to  calculate  the  development  of  the 
state  densities  in  an  expanding  flow.  All  that  is 
required  for  calculation  are  the  appropriate  rate 
coefficients  and  the  appropriate  thermodynamics 
for  the  molecular  internal  distributions  chosen. 
With  these  physical  properties,  the  conventional 
finite-rate  chemistry  gas-dynamic  code  can  be 
used  to  calculate  the  internal  distribution  of  states 
as  the  gas  flow  develops  along  the  nozzle  as  if 
each  state  were  a  separate  chemical  specie.  A 
separately  defined  vibrational  temperature  for  the 
description  of  the  internal  distribution  function  is 
unnecessary.  In  this  manner,  a  non-Boltzmann 
vibrational  distribution  will  manifest  itself  as  a  natu¬ 
ral  result  of  the  calculation  when  the  collision  fre¬ 
quency  becomes  too  low  to  maintain  a  particular 
state  in  collisional  equilibrium  with  the  rest  of  the 
ensemble.  This  is  the  approach  taken  here  and  the 
details  of  the  molecular  model  and  thermodynam¬ 
ics  used  for  this  calculation  are  described  below. 

Thermodynamics 

The  flow  code  used*’’*"'^^  requires  curve  fits  in 
the  NASA  form''^  for  the  species  nondimensional 
enthalpy,  H/RT,  specific  heat,  Cp/R,  and  entropy, 
S/R.  These  quantities  are  defined  for  an  ensemble 
gas  at  thermal  equilibrium  and  Include  the  contribu¬ 
tions  of  each  of  the  states  considered  here.  Since 
the  molecular  approach  here  is  to  treat  each  sepa¬ 
rate  collection  of  molecular  states  as  a  separate 
species,  it  was  necessary  to  generate  curve-fit 
coefficients  for  each  of  the  collections  considered. 


Str 

=  1.1512(5logioT  +  3logio*^)-‘>'‘>®4^ 

and  M  Is  the  molecular  weight,  T  is  the  temper¬ 
ature,  and  G(vo)  is  the  reference  zero  point  energy 
for  the  manifold.  The  Internal  partition  functions  are 
calculated  only  over  those  specific  vibrational 
states  comprising  the  particular  manifold: 


Qint  =  1  Qr,(T)gve 

''0 

Although  the  rotational  states  are  assumed  to 
be  in  thermal  equilibrium  at  the  gas  temperature, 
the  rotational  partition  function  was  calculated  for 
each  vibrational  state.  The  molecular  constants  for 
the  molecular  potential  functions  were  taken  from 
Ref.  15  for  the  diatomic  molecules  N2  and  O2  and 
Ref.  16  for  C02-  Since  the  rotational  constant,  Bv, 
and  the  centrifugal  constant,  Dq,  may  be  different 
for  each  of  the  different  stretching  and  bending 
modes,  they  were  calculated  for  each  level. 

The  thermodynamic  properties  for  each  collec¬ 
tion  of  states  were  calculated  in  100  K  increments 
from  200  to  3,000  K  and  the  results  least  square 
curve  fit  to  the  functional  form  of  the  NASA  coeffi¬ 
cients.  Two  sets  of  coefficients  were  generated, 
one  set  for  temperatures  below  1 ,000  K  and  the 
other  set  for  the  temperature  range  from  1 ,000  K  to 
3,000  K.  The  curve  fits  were  generated  over  both 
temperature  ranges  simultaneously  and  were  con¬ 
strained  to  be  continuous  and  smooth  at  1 ,000  K. 

RATE  COEFFICIENTS 


This  was  accomplished  with  the  specific  heat, 
enthalpy,  and  entropy  defined  from  the  partition 
functions  in  the  usual  way.^^ 

Cp  5  d  ^  2d//7QinA 

R  2  dn  dT  j 

H  _  5 

RT  2^^  RT  dT 


where 


S 

R 


R 


+  /'jQint 


+  T 


d//7Qin, 

dT 


The  molecular  model  for  CO2  required  for  this 
study  Included  significantly  higher  vibrational 
states  than  had  been  included  in  previous  work 
(e.g.,  Ref.  6),  and  a  different  distribution  of  individ¬ 
ual  microstates  within  each  discrete  collection  was 
used  for  the  computation.  Consequently,  it  was 
necessary  to  review  the  rate  coefficients  used  pre¬ 
viously  for  application  to  the  present  problem.  For¬ 
tunately,  considerable  experimental  work  was  per¬ 
formed  on  the  state-specific  relaxation  rates  of 
CO2  during  the  late  1960’s  and  early  1970’s  and 
significant  compilations  of  rates  exist.  The 
measurements  of  that  era  concentrated  on  the 
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states  in  the  V3  =  0  manifold,  ranging  from  the 
ground  state  to  (05®0),  and  there  has  been  littie 
work  in  this  area  since.  Accurate  modeling  of  the 
CO2  reiaxation  in  the  RDHWT  depends  greatly 
upon  energy  transfer  between  very  highly  excited 
vibrational  states,  including  those  much  higher 
than  (05®0)  and  other  V3  manifolds.  Thus,  it  was 
necessary  to  define  a  method  of  incorporating 
measured  energy  transfer  rates  where  availabie, 
and  generate  reliable  reaction  rates  for  unmea¬ 
sured  processes. 


identical  temperature  dependence  as  the  lower 
state.  This  crucial  assumption  then  alleviates  the 
need  for  detailed  first-principle  calculations  of 
velocity-dependent,  state-specific,  collision  cross 
sections  for  the  large  number  of  upper  CO2  states. 
The  variation  of  the  VT  rate  for  reactions  within  a 
given  subset  will  behave  as  the  variation  of  the 
transition  matrix  element,  which  is  a  function  of  the 
initial  set  of  vibrational  quantum  numbers.  Nine¬ 
teen  subsets  of  reactions  with  defined  vibrational 
quanta  changes  were  identified  in  this  analysis. 


Figure  5  shows  a  small  portion  of  the  CO2  term 
diagram  of  Fig.  3  and  includes  some  of  the  detailed 
states  of  the  V3  manifold.  Denoted  in  Fig.  5  are  the 
transitions  in  which  the  rate  was  either  measured 
(solid  line)  or  theoretically  calculated  (dashed  line). 
As  mentioned  above,  the  measured  rate  data  is 
complete  for  all  transitions  shown,  including  the 
ground  state  through  the  (05®0).  To  determine  the 
VT  rates  for  the  other  upper  states  required  here, 
an  extrapolation  procedure  was  developed  which 
was  based  upon  the  separability  of  the  dynamic 
collision  process  from  the  quantum  mechanical 
probability  for  transition  between  eigenstates.  The 
set  of  VT  reactions  can  be  broken  into  subsets,  in 
which  each  subset  is  composed  of  those  relaxation 
reactions  which  undergo  a  specified  change  in 
vibrational  quanta.  Since  the  collision  dynamics 
are  a  consequence  of  the  inelasticity  of  the  colli¬ 
sion,  each  reaction  of  a  given  subset  will  have 
identical  collisional  behavior  because  of  the  identi¬ 
cal  change  in  energy  and  angular  momentum. 
Hence,  VT  reactions  of  upper  states  can  be  scaled 
from  lower  states  within  the  same  subset  with  an 


The  calculation  of  the  transition  matrix  must 
consider  the  Fermi  resonance  of  the  vi  and  V2 
modes.  This  accidental  degeneracy  leads  to  a  mix¬ 
ing  of  states  of  similar  symmetry  and  the  resultant 
state  is  denoted  by  parentheses  enclosing  all 
degenerate  states.  Fpr  example,  the  degenerate 
states  (02°0, 10°0)  are  written,  to  first  order,  as^ 

(02°0, 10°0)  =  0.73  I  02°0>  +  0.68  I  10°0> 

and 

(02°0, 10°0)  =  -0.68  I  02°0>  +  0.73  I10°0> 

The  state  expansions  are  Important  in  the  cal¬ 
culation  of  the  transition  matrix  elements  between 
states  since  the  coupling  coefficients  lead  to  addi¬ 
tional  nonzero  transition  elements  due  to  the  cross 
terms  of  the  interaction  Hamiltonian.  For  the  upper 
states  of  the  V3  manifolds,  the  degree  of  the  Fermi 
degeneracy  increases,  such  as  for  the  (08^0, 16^0, 
24^0,  32^0)  state.  Thus,  relaxation  of  this  state  Into 
the  (07^0,  15^0,  23^0)  can  invoke  numerous  cross 
terms.  In  this  analysis,  the  eigenstates  were  com- 


I  Data 

I 

I 

^  Theory 


Fig.  5.  Existing  rate  constants  for  C02  VT  reactions. 
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puted  to  second  order  in  order  to  capture  the  third- 
order  interaction  amplitude.  This  degree  of  expan¬ 
sion  proved  especially  useful  when  incorporating 
relaxation  transitions  involving  the  v3  mode. 

As  an  example  of  the  rate  scaling  methodology, 
consider  two  subsets  of  relaxation  as  denoted  on 
Fig.  6,  denoted  as  Group  I  and  Group  II.  In  Group  I, 
the  set  of  relaxation  processes  is  characterized  by 
the  change  of  spectroscopic  quantum  numbers 
(0  n  "-2  0, 1  n-2  0)  (0  n-1  0, 1  n-3  0). 

Note  that  for  this  subset  of  relaxation  reactions, 
three  reactions  were  measured, 

(03^0, 11’0)^(02°0, 10°0), 

(0420, 1220)^(03^0, 11'’0), 

(OS^O,  1320)  ->  (0420,  122o), 

as  indicated  by  the  solid  lines  in  Fig.  6.  The  rates 
were  tabulated  in  the  standard  Landau  format. 
Irk  =  A  +  B/T  +CT''^2  jp  Rgf  ■\  s.  Figure  7  shows  the 


value  of  the  parameter  A  for  these  three  reactions. 
Also  shown  in  this  figure  is  the  extrapolation  of  the 
relaxation  rate  parameter  A  for  relaxation  of  higher 
vibrational  states  of  the  same  subset.  The  variation 
of  the  A  parameter  is  well  modeled  for  those  rates 
where  measured  data  exists,  remembering,  of 
course,  that  one  rate,  in  this  case  the  lowest  rate, 
was  used  to  scale  the  absolute  magnitude  of  the 
extrapolation. 

The  Group  I  rates  shown  in  Fig.  6  are  domi¬ 
nated  by  the  first-order  term  in  the  interaction 
amplitude,  i.e.,  the  change  of  only  one  vibrational 
quanta.  The  more  complicated  Group  II  transitions, 
which  involve  V3  mode  transitions,  require  second- 
and  third-order  terms  for  correct  physical  represen¬ 
tation,  i.e.,  v,  and  V2  transitions  as  well.  Since  the 
computation  of  the  third-order  terms  required  insig¬ 
nificant  additional  effort,  all  reaction  rates  were 
eventually  calculated  to  this  order,  regardless  of 
the  magnitude  of  the  third-order  term  relative  to  the 
lower-order  terms. 


Group  I 


0680- 


0770- 


0600- 


0650- 


OMO- 


0830“ 


0220- 


0110- 


Group  II 

--^i- 


0990  :(O07O,177D)  (0950.1750,2550^o980.1730,2^Q,a83C^^10.1710^^  10) 


-0561- 


0441- 


:(0750,1560)(0730,1530^)  (0710^0,2310;3110) 
1!o640,1440](0620, 14202^  (0600J^,22003000) 


0221- 


^^5;^j(0510.13102110) 

||^^^i^,12P0g000) 


^^^11- 


(0310,1110) 


(1000,0200) 


'  (0641,1441)  (0621,1421,2221)  (0601,1401, 
2201,3001) 


(0531,1331)  (0511,1311,2111) 


(0421,1221)  (0401,1201,2001) 


(0311,1111) 


(1001,0201) 


•0001 


Group  1-Data 


Group  I-Extrapolated 


Group  II-  Data 


Group  ll-Extrapolated 


Fig.  6.  Example  extrapolation  groups. 
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Fig.  7.  Data  and  extrapolated  rate  constants,  A,  for 
Group  1  reactions  of  Fig.  6. 

With  the  state-specific  relaxation  rates  extrapo¬ 
lated  from  the  measured  data  for  almost  all  the 
rates  of  interest,  some  method  of  reducing  the 
number  of  total  CO2  states  to  coincide  with  the  LDD 
modeling  was  necessary.  With  the  collecting  of  the 
individual  microstates  into  a  single  state  described 
above,  the  final  issue  is  to  determine  the  proper 
combination  of  state-specific  rates  that  describes 
the  relation  between  macrostate  relaxation. 

Consider  a  macrostate  N,  defined  by  a  certain 
combination  of  true  states.  The  number  density  of 
the  macrostate,  N  j,  is  then  the  sum  of  the  number 
densities  of  the  true  states,  N  j.  The  time  variation  of 
the  true  state  densities  may  be  given  as. 


macrostate  are  defined  by  thermal  equilibrium.  This 
is  reasonable  if  the  energy  differences  between  the 
states  are  small,  and  constitutes  the  driving  rule-of- 
thumb  for  macrostate  definition.  In  this  instance,  the 
macrostate  rate  can  then  be  written  as 

-M 


-E 

_ L 

kT 


kL„  = 


M^N 

U 


where  Qm(T)  denotes  the  partition  function  taken 
over  the  M  macrostate.  Thus,  the  net  relaxation  rate 
between  macrostates  can  be  written  as  a  properly 
averaged  sum  of  the  state-specific  relaxation  rates 
between  all  the  components  of  the  two  mac¬ 
rostates.  This  computation  was  carried  out  for  the 
19  VT  macrostate  relaxation  processes  used  in  the 
flow  computation. 


Results 

The  gas  dynamic  results  of  the  calculations  are 
shown  in  Figs.  8-10.  Figure  8  is  the  Mach  number, 
Fig.  9  is  the  static  temperature,  and  Fig.  10  is  the 
static  pressure.  Results  of  four  calculations  are 
included  on  each  figure: 

1 .  An  LTE  calculation  with  no  heat  addition 

2.  A  LTE  calculation  with  heat  addition 


3.  A  state  specific  relaxation  (SRR)  calculation 
with  no  heat  addition 

4.  An  SSR  calculation  with  heat  addition 


where  X  denotes  the  arbitrary  collision  partner  and 
Mj  denotes  states  from  macrostate  M  which  relax 
into  macrostate  N.  The  rate  of  change  of  the  density 
of  the  entire  macrostate  follows  as 

Nt  =  SNi 

]  i  ^  i  j  ' 

=  xfK' 

MN  '  N  V 

where  the  final  line  serves  as  a  definition.  To  allow 
for  the  macrostate  rate  solution,  an  assumption  is 
made  that  the  densities  of  the  true  states  within  a 


Additionally,  on  Fig.  9,  the  static  temperature  for 
a  calculation  in  which  the  10-kW  input  power  was 
added  to  the  plenum  chamber  rather  than  in  the 
downstream  supersonic  portion,  is  shown.  On  the 
scale  of  these  figures,  there  is  no  discernible  differ¬ 
ence  between  the  LTE  and  SSR  calculations  when 
there  is  no  energy  addition  in  the  supersonic 
stream.  Hence,  the  lines  corresponding  to  calcula¬ 
tions  1  and  2  above  are  indistinguishable.  The  LTE 
calculations  were  accomplished  with  the  same 
state  reactions  used  for  the  SSR  calculations, 
except  with  the  rate  coefficients  increased  to  ther¬ 
mal  rates.  That  is,  the  rate  coefficient  is  increased 
so  that  every  molecular  collision  results  in  a  transi¬ 
tion.  Included  also  on  each  of  the  figures  is  a  nor- 
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static  Pressure,  Atm  Static  Temperature,  K  Mach  Number 


Fig.  8.  Mach  number  axial  profiles. 


Fig.  9.  Static  temperature  axial  profiles. 


Fig.  10.  Static  pressure  axial  profiles. 


malized  energy  addition  profile  (Fig.  2) 
for  easy  reference. 

Clearly,  the  energy  addition  and  the 
vibrational  relaxation  has  an  effect  on 
the  flow.  Examining  Fig.  9,  it  is  seen  that 
adding  the  energy  to  the  plenum  cham¬ 
ber  results  in  an  increase  in  the  LTE  exit 
plane  temperature  of  about  25  K  above 
the  LTE  with  no  heat  addition.  There  are 
attendant  changes  in  Mach  number, 
velocity,  temperature,  etc.  For  the  calcu¬ 
lations  in  which  the  power  input  is  down¬ 
stream  of  the  throat,  all  calculated  prop¬ 
erties  are  coincident  upstream  of  the 
energy  addition.  Shortly  downstream  of 
the  point  at  which  the  energy  addition 
commences,  the  SSR  and  LTE  calcula¬ 
tions  begin  to  show  deviations  from  the 
baseline  curve.  For  the  region  of  energy 
addition,  the  LTE  calculation  shows  a 
higher  temperature  and  lower  Mach 
number  than  the  SSR  calculation,  as 
one  would  expect.  Since  the  LTE  calcu¬ 
lation  thermalizes  the  input  energy 
almost  immediately,  that  energy  is 
available  for  heating  of  the  flow.  In  the 
case  of  the  SSR,  the  input  energy  is 
maintained  in  the  internal  distribution 
and  is  not  released  to  the  gas  until  fur¬ 
ther  downstream.  Interestingly,  at  the 
exit  plane,  the  LTE  calculation  shows 
higher  Mach  number  and  lower  temper¬ 
ature  than  the  SSR,  indicating  that 
there  is  still  significant  energy  stored  in 
the  vibrational  distribution  of  the  mole¬ 
cules.  The  static  pressure  shows  similar 
characteristics,  with  the  LTE  calculation 
showing  the  greatest  increase  above 
the  baseline  calculation  and  the  SSR 
and  LTE  calculation  crossing  near  the 
exit  plane. 

Because  there  is  a  demonstrable  dif¬ 
ference  at  the  exit  plane  between  the 
SSR  and  LTE  calculations,  it  is  of  inter¬ 
est  to  examine  the  internal  energy  dis¬ 
tribution  determined  by  the  SSR  calcu¬ 
lation.  This  is  shown  in  the  form  of  Boltz¬ 
mann  plots  (Figs.  11  and  12)  of  the 
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Fig.  1 1 .  Exit  plane  Boltzmann  plot  of  CO2  Internal  state 
densities. 


Fig.  12.  Exit  plane  Boltzmann  plot  of  selected  CO2,  N2, 
and  O2  densities. 

excited  species  at  the  exit  plane.  The  ordinate  is 
the  logarithm  (base  10)  of  the  number  density  of 
the  specific  vibrational  state  relative  to  the  statisti¬ 
cal  weight  and  the  ground-state  density,  and  the 
abscissa  is  the  internal  energy.  On  a  plot  of  this 
nature,  those  states  that  are  in  thermal  equilibrium 
with  each  other  lie  on  the  same  straight  line  at  a 
slope  of  (-1/kT)  where  k  is  Boltzmann’s  constant 
and  T  is  the  temperature  of  the  distribution. 

Figure  1 1  shows  the  Boltzmann  plot  at  the  exit 
plane  of  the  CO2  internal  states  only  and  the  NLTE 
nature  of  the  distribution  is  evident.  Each  state 
included  for  CO2  is  shown  except  for  the  artificial 
upper  manifolds  of  states.  Examining  the  detailed 
states  plotted,  clearly  there  is  no  single  straight 
line  that  can  be  used  to  characterize  the  CO2 
vibrational  distribution.  There  are  some  apparent 
relationships  between  collections  of  states  that  are 
worth  examining  and  there  are  three  lines  drawn 


on  the  figure  to  illustrate  these  relationships.  It 
is  emphasized  that  the  only  temperature  used 
throughout  the  calculations  was  the  gas  transla¬ 
tional  temperature,  and  it  is  shown  and  identi¬ 
fied  on  the  plot  as  Tgas  =  198K.  The  states  with 
specific  V3  appear  to  fall  on  the  straight  line  at 
some  other  temperature.  This  is  illustrated  on 
the  figure  by  a  line  drawn  through  the  states 
with  V3  =  0.  This  equivalent  temperature  is  at 
250  K,  and  is  significantly  different  than  the  gas 
temperature.  Although  nearly  at  the  same 
slope,  the  distributions  for  the  series  of  states 
V3  =  1 , 2  appear  to  be  at  a  slightly  different  tem¬ 
peratures  compared  to  each  other  and  to  the  V3 
=  0  series.  This  difference  in  slope  may  easily 
be  a  manifestation  of  the  effect  of  the  assump¬ 
tions  about  the  upper  states  and  their  coupling 
to  the  next  lower  state.  If  enough  states  were 
included  in  each  manifold,  the  temperature 
characteristic  of  the  distribution  of  states  within 
each  manifold  might  easily  be  the  same.  The 
lowest  states  in  each  of  the  V3  manifolds 
appear  to  fail  on  the  same  straight  line  also, 
albeit  there  are  only  three  states  to  determine 
this.  The  apparent  equivalent  temperature  for 
this  relationship  is  1 ,000  K.  Clearly,  a  descrip¬ 
tion  of  the  CO2  internal  distribution  will  require 
accounting  for  at  least  three  temperatures. 

Figure  12  shows  an  exit  plane  Boltzmann 
plot  Including  the  N2  and  O2  vibrational  states  but 
including  only  the  CO2  (000),  (001),  and  (002) 
states.  The  lines  at  the  temperatures  250K  and 
1 ,000  K  are  shown  for  easy  reference.  The  strong 
coupling  between  the  states  of  N2,  and  the  V3 
modes  of  CO2  is  evident.  The  N2  vibrational  states 
nearly  overlay  the  (001)  and  (002)  states  of  CO2 
and  appear  to  exhibit  the  same  distribution,  as  evi¬ 
denced  by  the  straight  line  at  1,000  K  drawn 
through  those  densities.  The  O2  vibrational  states, 
however,  exhibit  a  different  relationship  which  has 
a  lower  characteristic  temperature,  here  about  600 
K.  This  is  a  consequence  of  the  N2  -  O2  V-V  relax¬ 
ation  rate  and  the  V-T  relaxation  of  O2  being  at  dif¬ 
ferent  characteristic  rates  than  the  N2  -  CO2  V-V 
and  V-T  relaxation. 

Examination  of  Boltzmann  plots  at  different 
axial  positions  shows  essentially  the  same  results, 
with  the  V3  manifolds  of  CO2  exhibiting  a  tempera- 
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ture  above  the  gas  translational  temperature,  the 
N2  and  C02(00v3 )  exhibiting  another  temperature, 
and  the  O2  exhibiting  yet  another.  It  appears  that, 
for  accurate  description  of  the  energy  deposition 
and  the  relaxation  of  the  molecules,  at  least  four 
temperatures  will  be  necessary.  It  is  notable  that 
the  C02(02L1)  state,  which  absorbs  all  the  energy 
of  the  beam,  does  not  show  significant  departures 
from  the  distribution  with  its  neighboring  states  with 
the  same  V3,  even  in  the  axial  region  of  greatest 
power  addition.  It  should  be  noted  that  there  is  in 
fact  an  eievation  of  this  state  above  the  distribu¬ 
tion,  but  it  is  not  discernible  on  plots  of  this  scale. 
This  implies  that  the  thermalization  of  the  energy 
into  the  adjacent  states  is  in  fact  quite  rapid,  and 
the  (02L1)  state  never  shows  a  great  departure 
from  its  instantaneous  internal  distribution  which  is 
determined  by  the  relaxation  of  CO2  as  a  whole. 

It  is  instructive  to  examine  the  relative  mole 
fraction  of  several  of  the  states  important  in  the 
relaxation  process.  This  is  shown  in  Fig.  1 3,  which 
shows  the  axial  profiles  of  the  relative  mole  frac¬ 
tions  of  the  (000)  ground  state,  the  (02L1)  state 
which  is  pumped  by  the  laser,  and  the  (001)  state, 
the  lowest  state  in  the  manifold  including  the 
(02L1)  state.  The  molecular  mole  fractions  have 
been  normalized  by  their  maximum  value  so  that 
relative  comparisons  can  be  made.  Here  it  can  be 
seen  that  the  (02L1)  state  shows  a  very  rapid 
increase  in  mole  fraction  with  the  heating  profile 
increase.  At  the  same  time,  since  the  laser  is 
pumping  molecules  from  the  ground  (000)  state  to 
the  (02L1)  state,  the  (000)  state  shows  a  slight 
decrease.  The  relative  changes  in  the  two  curves 
is  significantly  different  because  there  are 
about  four  orders  of  magnitude  more  mol¬ 
ecules  in  the  ground  state  than  in  the  c  ^  q 
(02L1)  state.  The  (02L1)  state  continues  •■§ 
to  increase  in  mole  fraction  untii  about  25  £  0.8 

cm,  at  which  there  is  rapid  decay  consis-  ■§ 
tent  with  the  assumed  decrease  in  energy  ®  °'® 

deposition  from  the  iaser  beam.  The  (001)  J 
state,  which  is  connected  only  collisionally  | 
to  the  (02L1)  state  through  the  rate  coeffi-  J  0  2 
dents,  shows  an  increase  similar  to  the 
(02L1)  state,  albeit  iagging  slightly,  Down-  ® 
stream  of  the  heating  profiie  the  (001)  and 
(02L1)  states  appear  to  decay  exponen- 
tiaiiy,  reflecting  the  deactivation  of  the 


species  and  the  transfer  of  the  internal  energy  into 
the  translational  temperature.  The  (000)  state, 
however,  shows  an  increase  in  mole  fraction  as  all 
the  molecules  that  were  pumped  into  higher  states 
decay  back  down. 

It  is  of  interest  to  examine  the  energy  content  of 
the  vibrational  distribution  relative  to  the  energy 
input  to  the  gas  and  included  in  the  excited  distri¬ 
bution.  This  is  shown  in  Fig.  14.  Three  results 
curves  are  shown  on  Fig.  14  in  addition  to  the  rela¬ 
tive  energy  addition  profile.  Two  of  the  cunres  illus¬ 
trate  the  total  power  injected  into  the  (02L1)  state 
by  the  laser  and  the  power  carried  in  the  excited 
NLTE  internal  energy  states  in  excess  of  that  in  the 
LTE  distribution,  respectively.  These  curves  use 
the  right  ordinate.  The  other  curve  is  the  fractional 
portion  of  the  input  power  that  is  stored  in  the 
NLTE  distribution  beyond  the  LTE  content  of  the 
excited  states  and  uses  the  ieft  ordinate.  The 
power  in  the  excited  NLTE  distribution  is  deter¬ 
mined  by  simply  adding  up  the  number  of  mole¬ 
cules  in  each  state  multiplied  by  the  energy  of  that 
state  and  dividing  by  the  instantaneous  velocity 
and  area.  The  power  in  the  LTE  distribution  is 
arrived  at  the  same  way,  except  that  the  population 
density  of  each  state  is  caiculated  from  the  Boltz¬ 
mann  distribution  at  the  instantaneous  gas  static 
temperature.  The  totai  power  injected  into  the 
(02L1)  state  is  simply  the  integral  of  the  energy 
addition  profile  shown  in  Fig.  2. 

Some  interesting  observations  arise  from  Fig. 
14.  First,  note  that  the  energy  content  in  the  NLTE 
distribution  rises  with  the  input  power  to  an  axial 
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Fig.  14.  Power  in  NLTE  distribution  calculated  from  SSR  model. 


position  of  about  25  cm,  at  which  point 
the  energy  deposited  into  the  gas  drops 
sharply.  At  that  point  collisional  pro¬ 
cesses  dominate  the  radiative  excitation 
process  and  there  is  a  decrease  in  the 
power  content  in  the  NLTE  distribution. 

This  is  consistent  with  physical  intuition, 
realizing  that  the  process  of  energy 
redistribution  to  the  other  modes  of  the 
gas  require  finite  times  to  accomplish 
whereas  there  is  an  instantaneous  per¬ 
turbation  in  the  excited  state  distribution 
because  of  the  laser  absorption.  The 
plot  of  the  excess  power  in  the  NLTE 
distribution  relative  to  the  total  power 
input  shows  an  almost  instantaneous 
increase  to  near  90  percent  at  an  axial  position  of 
about  10  cm,  the  point  at  which  energy  addition 
commences  in  the  flowing  gas.  The  relative  power 
then  starts  to  decay  slightly  to  an  axial  position  of 
about  15  cm  at  which  it  again  increases  to  near  87 
percent.  The  1 5-cm  axial  position  marks  a  point  at 
which  the  rate  of  energy  deposition  into  the  (02L1) 
state  increases  significantly.  Near  10  cm,  the 
instantaneous  population  of  the  (02L1)  state 
becomes  relatively  large  enough  that  the  collisional 
depopulating  processes  become  competitive  with 
the  rate  of  the  populating  laser  input.  Between  1 0 
and  15  cm,  since  the  rate  of  energy  input  is  con¬ 
stant,  the  collisional  processes  come  to  dominate 
the  NLTE  distribution  and  the  relative  power  in  the 
NLTE  distribution  decreases.  When  the  rate  of 
laser  power  input  increases  at  15  cm,  the  radiative 
process  again  dominates  the  collisional  for  a  short 
time,  and  the  power  in  the  NLTE  distribution  again 
increases.  From  this  point,  since  the  rate  of  power 
input  doesn’t  increase  again,  the  relative  power  in 
the  NLTE  distribution  continues  to  decrease.  At  an 
axial  position  of  about  25  cm,  the  rate  of  laser 
power  into  the  gas  has  a  sharp  decrease  and  the 
decay  characteristic  of  the  NLTE  distribution 
increases.  Significantly,  for  this  nozzle,  there  is  still 
about  19  percent  (note  the  left-hand  scale)  of  the 
input  laser  power  that  is  still  contained  in  the  NLTE 
distribution  at  the  nozzle  exit.  A  longer  nozzle  will 
allow  still  more  of  the  energy  to  be  converted  to  the 
kinetic  modes  of  the  gas  and  allow  larger  Mach 
numbers  at  the  nozzle  exit. 


Summary 

State-specific  relaxation  (SSR)  calculations  for 
a  demonstration  experiment  intended  to  show  effi¬ 
cacy  of  the  radiatively  heated  wind  tunnel 
(RDHWT)  concept  were  accomplished.  The 
RDHWT  is  a  one  of  several  proposed  methods  to 
add  energy  into  the  supersonic  portion  of  a  nozzle 
flow  to  increase  the  test  chamber  conditions  (static 
temperature,  Mach  number,  etc.)  to  hypersonic 
flight  conditions.  Material  limitations  in  current  wind 
tunnel  technology  preclude  test  times  of  more  than 
a  few  milliseconds  duration  and  with  air  composi¬ 
tions  representative  of  flight.  Adding  energy  into 
the  supersonic  portion  provides  the  opportunity  to 
increase  the  energy  of  the  flow  while  still  staying 
within  material  limits.  The  key  to  successful  wind 
tunnel  design  will  be  the  ability  to  predict  and  con¬ 
trol  the  energy  input  and  the  expansion  to  achieve 
a  priori  test  conditions.  To  this  end,  a  demonstra¬ 
tion  experiment  has  been  proposed  to  test  the  abil¬ 
ity  to  add  the  energy  in  the  supersonic  flow  with 
predictable  results.^ 

The  demonstration  model  investigated  operates 
at  much  lower  pressures  and  laser  powers  than  the 
full-scale  facility,  but  suffices  to  demonstrate  the 
heating  of  supersonic  air  in  a  predictable  and  con¬ 
trollable  manner.  The  working  gas  for  the  demon¬ 
stration  calculations  were  air  seeded  to  a  40-per- 
cent  concentration  of  COo.  The  laser  used  was  a 
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10-kW  HF  laser  which  will  couple  to  the  CO2  in  the 
flow,  pumping  the  ground  state  of  CO2  to  the  (02^1) 
state.  Plenum  pressure  was  about  100  atm  and 
temperature  was  1,150  K. 

The  SSR  model  treats  each  vibrational  level  of 
the  gas  molecules  as  a  separate  species  and  does 
not  depend  upon  any  of  the  usual  assumptions 
regarding  the  distribution  of  the  energy  internal  to 
the  molecular  structure.  With  this  approach,  the 
input  laser  energy  can  pump  a  specific  vibrational 
state  (as  it  will  in  an  experiment)  and  the  reiaxation 
from  that  specific  state  through  the  rest  of  the  dis¬ 
tribution  to  the  translational  modes  of  the  gas  can 
be  examined. 

Results  of  the  calculations  show  that,  at  the 
conditions  used,  most  of  the  deposited  energy  is 
thermalized  into  the  flow.  The  relaxation  has  an 
effect,  in  that  compared  to  a  thermal  equilibrium 
caiculation,  there  was  still  significant  energy  bound 
in  the  excited  NLTE  internal  distribution  at  the  exit 
plane.  Longer  nozzles  and  slower  expansions  wiil 
ailow  greater  thermaiization  of  the  input  power.  The 
relaxing  distribution  of  internai  vibrational  states  of 
the  flow  molecules  exhibited  definite  non-Boltz- 
mann  distributions,  although  states  with  certain 
excitation  characteristics  did  exhibit  thermal-like 
characteristics,  albeit  at  temperatures  other  than 
the  gas  translational  temperature.  For  CO2  the 
(OOV3)  states  appear  to  exhibit  one  temperature 
distribution  while  the  states  with  constant  V3  exhibit 
a  different  temperature  distribution.  The  CO2  (OOV3) 
and  the  N2(v)  states  exhibit  their  weil-known  strong 
coupling  and  manifest  the  same  thermal-iike  reia- 
tionship.  The  02(v)  states  also  appear  to  exhibit 
thermal-like  characteristics,  but  at  yet  another  tem¬ 
perature.  Correct  accounting  of  the  energy  content 
of  the  vibrationally  reiaxing  flow  will  have  to  include 
these  different  temperature-like  characteristics, 
since  significant  energy  can  be  stored  in  each  of 
the  different  modes,  and  the  details  of  the  release 
of  the  stored  energy  to  the  flow  will  affect  the  calcu¬ 
lation  of  the  exit  plane  properties. 
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